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C^ ■ Abstract 

D , The Fortuin-Kasteleyn mapping between the Ising model and the site-bond correlated perco- 

'■^ lation model is shown to be only one of an infinite class of exact mappings. These new cluster 

representations are a result of "renormalized" percolation rules correlated to entire blocks of 
spins. For example these rules allow for percolation on "virtual" bonds between spins not 

5^ I present in the underlying Hamiltonian. As a consequence we can define new random cluster 

theories each with its own Monte Carlo cluster dynamics that exactly reproduce the Ising 
model. By tuning parameters on the critical percolation surface, it is demonstrated numer- 
ically that cluster algorithms can be formulated for the 2-d and 3-d Ising model that have 
smaller autocorrelations than the original Swendsen-Wang algorithm. 



1. Introduction 

Our understanding of the Ising model has been greatly enhanced by the study of Coniglio- 
Klein [|l| critical percolation clusters and the classic work of Fortuin-Kasteleyn|2|, |3[ on the 
exact mapping of the Ising model into the random cluster model for site-bond correlated 
percolation. With the work of Swendsen and Wang Q this cluster representation has also 
lead to accelerated Monte Carlo dynamics with very small autocorrelation times ||5|, |6|, Q- 

In all of these endeavors the underlying clusters representation appeared to be unique — the 
probability of bond percolation, the only free parameter, being fixed by detailed balance or 
the requirement of scaling in the critical region. However, a central mystery of this beautiful 
formalism remains concerning the incompletely decoupling of clusters as manifest in divergent 
autocorrelations at the phase transition. It appears to us that a "perfect" representation of 
normal modes via clusters (analogous to fourier modes for a pure Gaussian theory) should in 
principle exhibit no critical slowing down (z = 0). 

Here we show that the underlying Fortuin-Kasteleyn map is not unique. 

By introducing longer range correlated percolation, we show that there are additional free 
parameters (in the null vector space transverse to the critical percolation surface) which are 
not fixed by the requirement of an exact mapping between the Ising and the random cluster 
partition functions. The percolation process is extended to longer length scales without leaving 
the critical surface on which the pair connectedness function coincides identically with the 
Ising correlation function. The new freedom to define "renormalized" critical clusters may help 
to understand the residual cluster interactions and lead to even more efficient Monte Carlo 
methods. By removing the rigid constraints of earlier methods, we also hope to formulate 
new cluster methods for more complex systems such as spin glasses [§]. Even generalizations 
to continuous molecular systems now appear less daunting. 

In the present article our new cluster methods will be explained exclusively in the context 
of the Ising model, although it should be obvious that they apply equally to all the recent 
extensions of the Swendsen- Wang algorithm such as those made by Brower and Tamayo||5| for 
the Landau-Ginzburg theory, by Wolff [^J for spin models and others [^, 0, |^- While shorter 
autocorrelations evidently can result, in all probability the dynamical universality class and 
the critical exponent z will not be changed unless we go to recursive methods that introduce 
virtual bonds on all length scales. Such extensions of our methods will be considered in future 
publications. 



The organization of this article is as follows. Section 2 presents the block percolation scheme 
through two simple examples, with illustrations of possible extensions. Section 3 defines 
the resulting cluster dynamics and computes autocorrelations and cluster distributions. The 
renormalized clusters are shown numerically to lead to reduced autocorrelations relative to 
the Swendsen-Wang method. Appendix A presents the general construction of our block 
correlated decomposition and the proof of equilibrium for the resultant cluster dynamics. 
Although the general argument is important both for its rigor and for the extraordinary 
range of alternative cluster methods it implies, the specific examples are probably easier to 
understand. 

2. Blocking Schemes for the 2-d and 3-d Ising Model 

Let us begin by studying in detail the simplest extension of the Fortuin-Kasteleyn (F-K) 
mechanism from blocks of 2 spins to blocks of 3 spins. This modest generalization already 
encounters several qualitatively new features. 

The conventional F-K mapping has two steps. The first step decomposes the Ising Hamiltonian 
into a sum over bonds (or 2 spin blocks) so that the Ising probability distribution is seen as 
a product of factors. 

plsing^^,-^ = - Y[ e-^(l-"'"^) (1) 

The second step decomposes each factor into two terms, 

g-/3(i-sis,) = Y^p^ A„(si, S2) = e-2/3 + (1 - e-''')6s„s, (2) 

n=0,l 

with the probabilities, 

Po = e-2^ Pl = l-e-2^ (3) 

uniquely determined by the requirement that probabilities for aligned vs unaligned spins be 
correct. The deletion (n=0) and percolation (n=l) terms have been designated by Ao(si, S2) = 
1 and Ai{si,S2) = Ssi,s2 respectively to accommodate a more general expansion for larger 
blocks. 

2.1 Three Spin Blocks 

The general expansion for a block of spins B takes the form, 

e-/3VB{..)= ^ p„A„(s,). (4) 



Each of the n terms is a unique product of delta functions representing one of the many possible 
ways to break it up into clusters. The only requirement is that the entire Hamiltonian be 
given by the sum over the set of blocks and that the expansion of the block into delta functions 
be complete. (See Appendix A for the proof and an even larger class of options.) 

For example the 3-spin block for a nearest neighbor Ising lattice (See Figure 1 and 2), B = 
(si, 52,53), can be split into 5 terms. 

However this most general 3-spin expansionQis no longer unique. Since there are 5 indepen- 
dent p's and only 4 independent equations, for (51,52,53) = (-|-,+,-|-), (—,+,+), (+,—,-!-), 
(-|-,-|-,— ), the linear equations must have a non-trivial null space. The general solution, 
Pi = PI + pnuii ^ jg ^j-^g sum of an original inhomogeneous F-K solution pf and one null 
vector, pf"" = (-1,1,1,1,-2). 

p^ = e-'^{l-e-'^) + C 

P2 = C 

p, = (I - e-'^)e-'^ + C 

p, = (1- 6-2/3)2 _2C, (g) 

This null vector expresses the simple fact that at least two Ising spins must be in the same 
state. 

{l-6s,,s,){l-Ss,,ss){l-Ss,,s,) = (7) 

No new solutions are present for the q state Potts model for q > 2 without going to larger 
blocks. The parameter C has been chosen to represent a virtual bond connecting 5i to S3 
so that at (^ = this solution reverts to the standard Fortuin-Kasteleyn decomposition. As 
the parameter ^ is increased the "real" bonds are reduced and the longer range "virtual" 
bond increased moving along the critical surface of the mapping. The remarkable simplicity 
of percolation relative to most critical phenomena is seen in the ability to find the critical 
surface analytically as a function of the critical Ising temperature. 

For Monte Carlo cluster algorithms, it will turn out that the only requirement is to make 
all the probabilities non negative, which in the neighborhood of the critical temperature 
restricts C to < C < e~^^''. For the 2-d Ising at the critical temperature, the limiting case. 



^Note that the identity, Ssi,s2 Ss2,s3 = <5si,s2 Ssi,s3 Ss2,s3, imphes that the last term is actuaUy a symmetric 
cluster of all three spins. 



C = e~^'= = (1 — v2)^ has the remarkable property that both the no percolation term pQ, and 
the fully percolated term p4, are exactly zero. Consequently any 3 spins must have precisely 
one percolated bond with 100% probability. 

To illustrate a little the range of options we have considered a variety of ways of partitioning 
the 2-d lattice into 3 spin (or two bond) blocks (See Figure 1). It is instructive to consider 
the limiting case (" = e~^^ for lattice la with the three spin blocks centered on all the "black 
sites". If we start with a checkerboard and use our limiting case again, (" = e~^^'=, the lattice 
will immediately percolate into two sublattices with infinite clusters on the "red" sites and 
single clusters on the "black" sites. The study of such extreme non-equilibrium configurations 
shows how different our clusters are relative to Coniglio-Klein site-bond correlated percolation 
clusters. 

It is easy to extend this example to the most general three spin blocks, 

e-/3l2(l-.1.2)-ft3(l-S2S3)-Al(l-.3.l) =p^ +p^ 5^^^^.^+p^ <5,„,3 +P3 5,,,,, + ^4 6,,,,, 6,,,,, (8) 

appropriate for continuous spin models on cubic or triangular lattices. Again the general 
solution can be expressed as pi = pf^ +p""", where the inhomogeneous solution is given by 
the F-K prescription applied to each bond separately and the null vector is unchanged. If 
we consider a red-black decomposition of the triangular lattice {(iij = (3) , the limiting case , 
C = e~^^'' = v3 at the critical point again has the remarkable property that one and only 
one bond is percolated with 100% probability. 

In the remainder of the paper we will consider in greater detail the resultant random cluster 
model and its Monte Carlo dynamics for the 3-spin block and plaquette blocking schemes. 
However as the reader will readily note, we have opened a Pandora box of alternative schemes. 
As the blocks are increased in size, there are more and more free parameters and neither the 
blocks nor the choice of these parameters have to be chosen to respect the symmetries of 
the underlying lattice. Consequently the exact constraint of mapping the Ising model into a 
random cluster scheme is far less restrictive than we might expected. 

2.2 Plaquette Blocking for the Ising Model. 

As a second additional example we briefly describe a scheme corresponding to 4 spin blocks 
or plaquettes. The energy can be expanded either, as a sum over all plaquettes dividing the 
interaction f3 for each bond by the number of ways it is shared (2 in 2-d or 4 in 3-d), or by 
breaking the lattice up into a sum over over red (or black) plaquettes on the lattice. As in 
the 3-spin block case the Boltzmann factor for each plaquette is decomposed in a weighted 



sum over products of delta functions A„ that constrain different subsets of spins inside the 
plaquette. 



P2 v''s3,siOsi,S2 ~r ^si,S2^S2,S4 + '^82,34054,33 + ^S4,s,j'Js3,si) + 

P3 \^Sl,S3"S2,S4 ~T~ "81,82^33,34) y Pa ^8\, 82^82,84^84, 8-),^8'3, 81 \ 

Pb ^8l,S4^82,SZ +P6 {^81,84 + ^82,83), (9) 

where si,S2-,s^ and S4 are the four spins that define the plaquette (see figure 3). Notice that 
the probabilities ps and pQ represent virtual percolation bonds across the diagonals that freeze 
spins which do not interact in the Ising Hamiltonian. 

This expansion is not the most general however, since we have imposed rotational symmetry 
in the spins (si, S2, S3, 54). For example the use of three spin blocks as shown in figure Id offers 
a legitimate example of a four spin expansion as well, which violates rotational symmetries on 
the plaquette. The most general four spin expansion has 15 parameters and only 7 constraints. 
The number of solutions grows faster than exponential as the size of the block increases — 
for a 5 spin block there are already 52 parameters but only 16 constraints. Our 5 spin star 
blocking (fig la) could be expanded with 36 adjustable parameters ! 

Imposing rotational symmetry for the Ising model, there are four independent linear equations 
for (si,S2,S3,'54) = (+,+,+,+),(-,+,+,+),(-,-,+,+),(-,+,-,+) and their solution is 
expressed again as the sum of the F-K inhomogeneous solution plus the null space — this 
time with three free parameters (Ci, C2; Cs) for three null vectors. 
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e-^P - 4C2 
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(10) 

The parameters (^ s are constrained by the condition that all probabilities must be positive 
numbers in the range (0, 1). These three parameters give us additional flexibility to tune 
the percolation process in different ways. For example, one can increase or decrease the 
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probability of diagonal bonds with respect to horizontal or vertical bonds and, as we will see 
in section 5, in this way change cluster properties and decorrelation times. The plaquette 
example can also be applied to the 3 state Potts model with the result that the rotational 
symmetric solutions has only one null vector, pT^"'' = (1,-1,2,1,— 6,1,— 1), expressing the 
fact that in a 3 state Potts model the four spins of a plaquette can all reside in different states. 
In general it appears likely that for a q-state Potts model one must resort to blocks bigger 
than q spins to find a non-trivial null space. 

Finally we would like to emphasize that rather elegant blocking schemes exist whereby bonds 
are shared between adjacent blocks with the energy divided between the two blocks. Even 
for the F-K map this is allowed due to the trivial identity on each shared bond that ((1 — 
p)^si,s2 +P) ^ ((1 ~P')^si,s2 +P') = (1 ~ PP')^si,s2 +pp' ■ For example in our blocking schemes, 
we can restore translational symmetry for the 2-d triangular and square lattices simply by 
the replacement /? — > |/3 throughout eq. ^ and eq. |l^ respectively. Each bond is split into 
two equal pieces. For the Plaquette percolation method this has the appealing feature that 
all the symmetries of the original theory are preserved with next-to-nearest neighbor virtual 
bond percolation allowed over the entire lattice. We feel that the more symmetric forms are 
more likely to express the proper scaling of critical clusters and thus lead to more efficient 
algorithms. 

3. Cluster Dynamics. 

In order to define a dynamics based on the percolation mappings described in the previous 
section one has to introduce for every blocking scheme, a set of auxiliary percolation variables 
ub = 0, 1, ... on each block. Then the Monte Carlo cluster dynamics is chosen to be a 
Markov process for the joint distribution p^"*" {si^ns), in which one alternatively chooses 
the percolation variables ub at fixed spins Sj and the spins for each cluster at fixed ub (see 
Appendix A). 

We outline our algorithm in the following steps. 

1. Choose a local blocking scheme in such way that all spin-spin interactions are accounted 
for in the ensemble of blocks. Expand each block into a fixed percolation representation 
to define the joint distribution, p^°*" (si,nB)- 

2. For each block obtain its spin configuration (sj) and choose a value for the "percolation" 



variable n = ub with normalized probability: 

PnA„(Si) 

"'"'"' = E„P„A„(..) <"' 

3. Label the connected components for the percolation graph using a cluster finding algo- 
rithm such as breadth-first search|lC], the Hoshen-Kopelman algorithm[f 1|, or a parallel 



algorithm [12|. 

4. Flip all of the spins in each cluster coherently with 50 % probability and return to step 
2. 

Basically this is very similar to the standard Swendsen-Wang algorithm except that the per- 
colation step involves probabilities correlated to entire blocks of spins. The variations on 
Swendsen-Wang suggested by Wolffj^] can be adapted to these clusters as well. There is a 
slight increase in complexity in the percolation step compared with Swendsen-Wang but the 
algorithm is still very straightforward. The values for the C,' s parameters can be taken to have 
different values for different block in the system. In our simulations we choose all of them to 
be fixed at a uniform value and then generate a look-up table of probabilities to be employed 
in the percolation process. The 3-spin and 4-spin block dynamics at Ci = are identical to 
the original Swendsen-Wang algorithm. 

3.1 Numerical Results 

Due to their simplicity we have implemented the 3-spin and plaquette blocking schemes as 
described in the previous section. In this section we analyze the autocorrelation times for 
equilibrium simulations and the dynamics of cluster distributions as the system relaxes to 
equilibrium. 

Autocorrelation times analysis 

All our simulations were carried out at the critical point of the infinite system. We studied 
the autocorrelation function for the squared magnetization as an average over Monte Carlo 
time (2.5 x 10^ steps). Figure 4 shows two autocorrelation functions for a system of size 
64 X 64 using the 3-spin block dynamics with the star blocking shown in figure la. The two 
curves correspond to the values of C, at the two extremes of the interval (0, e~'^^'') of allowed 
values. The two functions show exponential relaxation but the extreme case farthest from 
the Swendsen-Wang point is faster by a noticeable amount. Moreover, the data in figure 5 
appears to exhibit a monotonic decrease in autocorrelation times as the probability of virtual 



bonds increases. A similar picture appears for the 3-spin block dynamics in three dimensions 
where for C = e'^^", r^2 is almost twice as fast as Swendsen-Wang. 

For the red-black plaquette dynamics we explored different regions of the solution space. We 
have not done an exhaustive search but it is clear that the autocorrelation times decrease 
as Ci decreases and the other two parameters, C2 and ^3, increase. The fastest dynamics, 
the one with smallest autocorrelation times, we were able to find is the one corresponding 
to Ci = -0.056349, C2 = 0.00245, Cs = 0.0098 which is about 1.6 times faster than the 
Swendsen-Wang dynamics. Similar increases in speed are obtained for the symmetric form 
which expands both red and black plaquettes simultaneously using eq. 10 with f3 — > ^P- 

A general trend is seen in all these cases, as one might expect, increasing the probability of 
long range virtual bonds reduces the autocorrelation times. 

Figure 6 shows the scaling of autocorrelation times with system size for the 3-spin block, 
plaquettes and Swendsen-Wang dynamics in two dimensions. If indeed the scaling is the 
same for all these dynamics, then they will have the same dynamic exponent z as the original 
Swendsen-Wang and consequently would belong to the same dynamic universality class [^,[T^. 
A more comprehensive and systematic study of these dynamics is required to settle this 
question but all our numerical results appear to be consistent with universality. 

Cluster relaxation analysis 

In addition to the standard autocorrelation times analysis we studied the equilibrium cluster 



distributions and the dynamics of cluster relaxation (see Stauffer, Kertesz and Miranda |14|). 
Figure 7 shows the equilibrium cluster distributions for a 32 x 32 system using the plaquette 
blocking and Swendsen-Wang. Notice that the two cluster distributions are not identical but 
have very similar scaling behavior. The fastest plaquette dynamics having larger probabilities 
for virtual bonds favors larger clusters. It is clear that the cluster distributions follow a power 
law with identical or similar exponents. These exponents are related to the fractal dimension, 
df, of the clusters by the relation, n^^ ~ g-('^+d/df) -^[^i^ h^q implication that our extended 
percolation clusters have the same or similar fractal dimension as the Coniglio-Klein clusters. 

In order to study cluster relaxation, a large number of simulations were performed in which 
the system relaxes to equilibrium at the critical point from an initial state with all spins up. 
Cluster histograms were accumulated at every time step and a logarithmic binning used to 
reduce fluctuations. Bin n, contains contributions from clusters sizes in the interval (2*, 2*+^ — 



1). Clusters numbers are assumed to relax exponentially to equilibrium, 

n,(t) ~ (1 - e-*/^0 <^. (12) 

The scaling of the cluster decorrelation time r^ with cluster size defines a critical exponent r, 

Ts ~ S" ~ s'/'^f (13) 

where z is the standard critical slowing down exponent and df is the fractal dimension of the 
clusters {df = d — jS/v). For the 2-d Ising model, dj = d — f3/v = 1.875. 

We have performed this analysis for two values of C, in the 3-spin block dynamics. We used 
a system of size 256^ and sampled over 400 experiments. Figures 8 shows the relaxation of 
ng{t) as a function of time for the Swendsen-Wang and the Q = e~^^'' dynamics. It is clear 
that the C = e"^^'' dynamics produces a much more rapid relaxation process for the cluster 
numbers, perhaps by a factor of 5. This is surprising in view of the much smaller change for 
the equilibrium autocorrelation time (a factor of about 1.5), and gives an indication that the 
non-equilibrium dynamics may depend very sensitively on the strength of virtual bonds. In 
fact, the relaxation is so fast that it is hard to extract the scaling behavior by computing 
relaxation times for different cluster numbers. A more thorough study of the effects of these 
dynamics on cluster properties would help us to understand how to tune the C parameters to 
have more control over the cluster relaxation process. 

All the numerical experiments confirm the validity of these dynamics as equilibrium Monte 
Carlo processes and their potential to reduce autocorrelation times. 

4. Conclusions 

The main purpose of this article was to demonstrate the existence of a new class of exact 
percolation mappings not allowed in the original methods of Fortuin and Kasteleyn. To 
clarify the discussion we have studied in detail 3 spin and 4 spin blocking schemes applied to 
the nearest neighbor Ising model on a cubic lattice. However as we demonstrate in Appendix 
A, the freedom to invent new cluster algorithms within this framework is extraordinarily 
unencumbered. Not only is the Fortuin-Kasteleyn mapping not unique, one can find exact 
mappings which have peculiar properties such as not respecting the exact symmetries of the 
underlying model. We also note as illustrated by eq. |^ that these new mappings can not 
in general be extended analytically to the entire q state Potts model. The lower the value 
of q the larger the class of solutions. Now the problem is no longer the difficulty of finding 
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percolation mappings but the problem of choosing the most interesting or efficient straw in 
the haystack. 

It is difficult not to be reminded of the analogous freedom in blocking schemes for the real 
space Renormalization Group (RG). In spite of this analogy several differences remain. For 
example the real space Renormalization Group has the following properties not shared by our 
approach: (i) the Hamiltonian has fewer degrees of freedom so that only the longer length 
scales are described correctly, (ii) the blocking procedure gives renormalized parameters that 
are automatically tuned to respect the scale change and (iii) the blocking is applied recursively. 
With respect to the first feature concerning the use of a semi-group, perhaps our approach 
offers an improvement in that we have an exact isomorphism. However the second two features 
we would like to emulate more closely, in the belief that they might be the key to better cluster 
methods. For example, it appears likely to us that the optimal choices for our Q parameters 
(ie those with shortest autocorrelations) can be "derived" from scaling arguments. We are also 
looking at recursive blocking schemes — in particular a class of recursive algorithms based 
on iterative block decomposition in the spirit of the Migdal-Kadanoff renormalization group, 
which may be able to reduce z below the value given by the Swendsen-Wang multi-cluster or 
Wolff single-cluster algorithms. Progress on these issues will be published in the future. 

It is also interesting to consider the possibilities that virtual bond percolation parameters 
may be able to tune percolation processes in cases in which the standard Swendsen-Wang 
percolation does not produce acceleration such as in spin-glasses. Work is also proceeding on 
this problem. In summary the use of virtual bond dynamics in its most general form may 
be able to open new perspectives for cluster methods both to understand critical domains or 
droplets and to design faster algorithms for Monte Carlo simulations in statistical mechanics 
and quantum field theory. 
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Figure Captions 

Figure 1. Examples of the variety of blocking schemes for 3-spins. The blocking used in 
our numerical experiments is shown in a). This blocking can be seen as a special case 
of a 5-spin "star" blocking. It is also interesting to notice that c) is an asymmetrical 
blocking and that d) produces plaquette-like structures. 

Figure 2. The expansion of the interaction between 3-spins in terms of percolation 
variables. 

Figure 3. Different percolation schemes for the 4-spin plaquette block mapping. 

Figure 4. Magnetization squared autocorrelation functions for the 3-spin block dynamics 
in two (a) and three (b) dimensions with different parameters C- The two values shown 
correspond to the extremes of the interval of allowed values. 

Figure 5. Variation of the magnetization squared autocorrelation time as a function of 
C for the 3-spin block dynamics in two dimensions. The value of t„2 decreases as C; the 
probability of virtual bonds, increases. 

Figure 6. Scaling of r„2 as a function of system size L for the 3-spin block, plaquettes 
and Swendsen-Wang dynamics in two dimensions. 

Figure 7. Equilibrium cluster distributions for Swendsen-Wang and plaquette dynamics 
(with Ci = -0.056349, C2 = 0.00245 and Cs = 0.0098) in two dimensions. The slopes are 
nearly identical but large clusters are favored by the plaquette dynamics. 

Figure 8. Relaxation of cluster numbers for Swendsen-Wang (a) and the 3-spin block 
(for C, = e"^^") (b) dynamics in two dimensions. 
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Appendix A: General Formulation of Block Percolation 

Our generalization of the Fortuin-Kasteleyn percolation map has a strong similarity to block- 
ing methods used for the real space Renormalization Group transformations. To bring out this 
similarity and to prove the correctness of the resulting dynamics, we follow the suggestion of 
Kandel et al[Q] by introducing effective block potentials, U'j^ {si,nB, /3), instead of the delta 
function constraints A„,^(sj) considered above. Also we note that the delta functions can be 
recovered as a limiting case. For example, in this approach the original F-K map is replaced 
by the following two terms, 

where C/^ (si,n,/3) = J„(l — S1S2) + log{pn)- If we adjust the two weights correctly, 

p^ = (e-2/3_e-2Ji)/(g-2Jo_g-2J,)^ 

this is in itself a (new) legitimate mapping of the Ising bond at temperature j3~^ into two 
randomly chosen terms with effective temperatures J^ . We can recover the standard F-K 
map by taking the limit of the effective inverse temperatures to zero ( Jq ^ 0) and to infinity 
(Ji -^ cxd) for the deletion and freezing (ie percolation) terms respectively (see eq.|2|). It is 
also interesting to view percolation methods in this way as the limiting extreme of a random 
mapping of a fixed temperature Ising system into a heterogeneous mixture of high and low 
temperature interactions. 

Our general blocking scheme has two steps: 

1. The Hamiltonian is expressed as a sum over blocks B of spins, 

H{si)=Y.VB{si). (A.3) 

B 

2. The corresponding factor in the partition function is decomposed into a sum of terms, 

nfl=0,l,...nmaa; — 1 

The products over all the resulting terms ( |A.^ ) forms a new joint probability distribution, 
P-'°*"*(sj,nB) so that our generalized map is 

{"s} 
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The Ising model is represented as a marginal distribution by summing over the generalized 
percolation variables, ns- Alternatively the sum over the Ising variables yields the generalized 
"random cluster" (RC) model, 

^''^(«b) = E^'""*(^^'^s). (A.6) 

Lastly by taking the limiting case for effective potentials to freeze or delete clusters of spins in 
each block, we can among others obtain the entire class of virtual bond (or extended) perco- 
lation models. The fundamental identity for this class of percolation models, that guarantees 
criticality of the random cluster model at the Ising fixed point, is 

(^S,s,) ^ Y. S^S,P''"'^{sk) = E l^,JP''^i^B) (A.7) 

{sk} {"b} 

where {siSj) is the Ising correlation function and the connectedness function, jij, is 1 if i and 
j are in the same random cluster and otherwise. From this identity it follows that the Ising 
spin-spin correlation length and the cluster connectedness length are equal and therefore that 
the critical scaling of both the Ising model and the corresponding random cluster model are 
identical. 

Equilibrium for cluster Monte Carlo 

To prove that all static quantities are given correctly by the Monte Carlo cluster algorithm 
in Section 3, one must show that the Markov process preserves the equilibrium distribution, 

plsing^^^^ = J2 W{s'i ^ Si) P^*^"3(si). (A.8) 

This equilibrium balance condition is weaker than the full detailed balance satisfied by many 
Monte Carlo algorithms, but it is sufficient. To prove it recall that in our cluster Monte 
Carlo algorithm the percolation variables ub in each block were first chosen according to the 
conditional probability, 

Pinsls.) = ^ = e-«'(-'-'/^)-/^^-(-)) (A.9) 

Z-^ub ^ 

and "clusters" of spins were subsequently chosen at fixed values of ub- Clearly this implies 
that the full transition matrix is a product of two conditional distributions, 

W{s[ ^Si)=J2 P-'°™*(s.|n's)P^""*(n'B|s.). (A.IO) 

Wb} 
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To complete the last step in demonstration one explicitly performs the sums in eq. A. 10 and 
eq. A.8| by applying Bayes theorem, 



Pia,b) = P{a\b)P{b) P(6)^^P(a,6), 

b 



(A.ll) 



for each of the marginal distributions in eq. A. 5 and eq. A.C. While there is nothing 



fundamentally new in this proof relative to that needed for the Swendsen-Wang algorithm it 
is necessary to see that it still goes through. Moreover this more general derivation illuminates 
a striking feature of cluster methods that the stochastic replacement of the full potential by 
a particular effective potential bears a very close resemblance to the Metropolis algorithm. 
Put into the context of blocks of spins, the similarity to the multigrid concept of "stochastic 
coarsening" is also transparent. 
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